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Abstract 

We investigate the causality and stability of relativistic dissipative fluid dynamics in the absence 
of conserved charges. We perform a linear stability analysis in the rest frame of the fluid and find 
that the equations of relativistic dissipative fluid dynamics are always stable. We then perform 
a linear stability analysis in a Lorentz-boosted frame. Provided that the ratio of the relaxation 
time for the shear stress tensor, r n , to the sound attenuation length, T s = 4n/3(e + P), fulfills 
a certain asymptotic causality condition, the equations of motion give rise to stable solutions. 
Although the group velocity associated with perturbations may exceed the velocity of light in a 
certain finite range of wavenumbers, we demonstrate that this does not violate causality, as long as 
the asymptotic causality condition is fulfilled. Finally, we compute the characteristic velocities and 
show that they remain below the velocity of light if the ratio r^/Tg fulfills the asymptotic causality 
condition. 
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I. INTRODUCTION 



Data from the Relativistic Heavy- Ion Collider (RHIC) on the collective flow of matter in 
nucleus- nucleus collisions have delivered a surprising result: the elliptic flow coefficient i>2 
is sufficiently large jl-4| to be compatible with calculations performed in the framework of 
ideal fluid dynamics Q|. This has given rise to the notion that "RHIC physicists serve up 
the perfect liquid" js-8]. 

Of course, no real liquid can have zero viscosity: for all weakly coupled theories, i.e., 
theories with well-defined quasi-particles, in the dilute limit there is a lower bound which 
one can derive from the uncertainty principle j^]: the ratio of shear viscosity to entropy 
density rj/s > 1/12. For certain strongly coupled theories without quasiparticles, there is 
also a lower bound which can be obtained from the AdS / CFT conjecture 10| , rj/s > 1/ (47r) , 
i.e., surprisingly close to the bound for dilute, weakly coupled systems. 

In order to see whether the shear viscosity of the hot and dense matter created in nuclear 
collisions at RHIC is close to the lower bound, one has to perform calculations in the 
framework of relativistic dissipative fluid dynamics. This program has only been recently 



initiated, but has already led to an enormous activity in the literature |12431|. 

Fluid dynamics is an effective theory for the long-wavelength, small-frequency modes of 
a given theory. In order to see this, let us introduce three length scales: (a) a microscopic 
length scale, ^ m i cro - In all theories, at sufficiently large temperatures this length scale can be 
defined as the thermal wavelength A t h ~ 1/T. In weakly coupled theories with well-defined 
quasi-particles, this can be interpreted as the interparticle distance, (b) A mesoscopic length 
scale, £ meS o- In weakly coupled theories and in the dilute limit, this can be identified with 
the mean-free path of particles between collisions. In strongly coupled theories, such a scale 
is not known and should be identified with £ m i CT0 . (c) A macroscopic length scale, £ maC ro- 
This is the scale over which the conserved densities (e.g. the charge density, n, or the energy 
density, e) of the theory vary. Thus, ^Lro ~ \de\/e, i.e., ^ m a Cro * s proportional to the 
gradients of the conserved quantities. 

We now define the quantity K = ^ m esa/^macro- For dilute systems, this quantity is identical 
to the so-called Knudsen number. If K is sufficiently small, fluid dynamics as an effective 
theory can be derived in a controlled way as a power series in terms K. Since K 
this series expansion is equivalent to a gradient expansion. 



^ macro' 
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To zeroth order in K, one obtains the equations of ideal fluid dynamics. To first order in 
K, one obtains the Navier-Stokes (NS) equations. So-called second-order theories contain 
terms of second order in K. Examples for the latter are the Burnett equations 32J, the 
Israel-Stewart equations for relativistic dissipative fluid dynamics [33} , the memory function 



theory 



26 



29|, extended thermodynamics 29|, |34j, and others [35|. The main difference 



between first and second-order theories is the velocity of signal propagation. The relativistic 
NS equations allow for infinite signal propagation speeds and are therefore acausal. On the 
other hand, all second-order theories are considered to be causal in the sense that all signal 
velocities are smaller than the speed of light, provided that the parameters of the theory are 
suitably chosen. 

The stability and causality of fluid-dynamical theories are usually studied around a hy- 
drostatic state (i.e., for vanishing macroscopic flow velocity) which is in thermodynamical 
equilibrium. However, if a theory is stable around a hydrostatic state, it does not necessarily 
imply that it is stable in a state of nonzero flow velocity. Following this idea, the stability and 
causality of first and second-order fluid dynamics for a state with nonzero background flow 



for the case of nonzero 
There it was found 



velocity (mathematically realized by a Lorentz boost) were studiec 
bulk viscosity, but for vanishing shear stress and heat flow in Ref. 
that causality and stability are intimately related: for all parameters considered, the theory 
becomes unstable if and only if there is a mode which propagates faster than the speed of 
light. 

In this paper, we extend this analysis to the case of nonvanishing shear viscosity in 
second-order theories of relativistic dissipative fluid dynamics. A similar analysis for a 



hydrostatic background has already been done by Hiscock, Lindblom, and Olson 



38 



but they discussed exclusively the low- and high-wavenumber limits [39|. As we shall show 
in this paper, their analysis missed a divergence of the group velocity of a shear mode at 
intermediate wavenumbers. This anomalous behavior is generic, i.e., it cannot be removed by 
tuning the parameters of the theory, e.g., the relaxation time for the shear stress tensor, t w , 
and the shear viscosity, 77. However, if the ratio T n /T s , where T s = 2(D — 2)rj/[(D — l)(e + P)] 
is the sound attenuation length in D space-time dimensions, is chosen such that the large- 
momentum limit of the group velocity associated with the perturbation remains below the 
velocity of light (the so-called asymptotic causality condition), one can ensure that the 
divergence is restricted to a finite range of momenta. It will be demonstrated that in this 
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case, the causality of the theory is not compromised. On the other hand, second-order fluid 
dynamics is always stable in the rest frame of the fluid, even if we use a parameter set which 
violates the asymptotic causality condition. 

We also study the causality and stability for a state with nonzero background flow velocity, 
i.e., in a Lorentz-boosted frame. We find that the divergence of the group velocity is removed. 
However, depending on the boost velocity the group velocity of either the shear or the sound 
mode may still exceed the speed of light in a certain range of wavenumbers. Nevertheless, 
provided that the ratio T n /T s fulfills the asymptotic causality condition, we can show that 
the equations are stable. In contrast to the analysis in the rest frame, however, they become 
unstable if the asymptotic causality condition is violated. We shall demonstrate that if 
the asymptotic causality condition is fulfilled, the causality of the theory as a whole is not 
compromised. In this sense, causality and stability are intimately related. 

So far, the discussion was limited to the fluid-dynamical equations in the linear approxi- 
mation. Therefore, we expect the results to be valid for all versions of second-order theories 
presently discussed in the literature, since they differ only by nonlinear terms. We also 
compute the characteristic velocities for the so-called simplified IS equations 16[ without 
linearizing these equations. Our analysis strongly indicates that the characteristic velocities 
remain below the velocity of light if the ratio T n /T s is chosen such that the asymptotic 
causality condition is fulfilled. 

The asymptotic causality condition implies that, for a given T s ~ rj, r n must not be 
arbitrarily small. This explains why relativistic NS theory is acausal, because there t w — > 0, 
while t] is non-zero. It also implies that second-order theories are not per se causal; they can 
violate causality (and become unstable) if a too small value for r n is chosen. The statement 
that second-order theories automatically cure the shortcomings of NS theory is therefore not 
true. 

This paper is organized as follows. In Sec. HIl we discuss the causality and stability 
of the linearized second-order fluid-dynamical equations in the local rest frame. We also 
extend this analysis to nonzero bulk viscosity and show that the divergence of the group 
velocity still exists in this case. In Sec. IHIt this discussion is generalized to a Lorentz-boosted 
frame. We discuss Lorentz boosts both in and orthogonal to the direction of propagation 
of the perturbation. It will be demonstrated that superluminal group velocities will not 
compromise the causality of the theory as long as the asymptotic causality condition is 
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fulfilled. In Sec. IIVt we compute the characteristic velocities in the nonlinear case. A 
summary of our results concludes this work in Sec. [V] An Appendix contains details of 
our calculations in Sec. [TV] The metric tensor is g^ u = diag(+, — , — , — ); our units are 
h = c = k B = 1. 



II. STABILITY IN THE REST FRAME 



As mentioned in the Introduction, there are several a ppr oaches to formulate a second 



order theory of relativistic dissipative fluids 



26 



2N 



29 



33 



35J . These approaches differ only 



by nonlinear (second-order) terms. However, since we shall apply a linear stability analysis 
in the following, these differences vanish and all approaches lead to the same set of linearized 
fluid-dynamical equations. In this work, we do not consider any conserved charges and thus 
are left with energy-momentum conservation, 



d^ v = , 



(1) 



where 



T"" = e u»u v - {P + n) A"" + tt^ (2) 

is the energy-momentum tensor. Here, e and P are the energy density and the pressure, 
while it M , n, and are the fluid velocity, the bulk viscous pressure, and the shear stress 
tensor, respectively. We also introduced the projection operator 



(3) 



which projects onto the (D — l)-dimensional subspace orthogonal to the fluid velocity. We 



compute in the Landau frame 



11] . where there is no energy flow in the local rest frame. 



In second-order theories of relativistic dissipative fluid dynamics, the bulk viscous pressure 
and the shear stress tensor are determined from evolution equations. In D space-time 
dimensions (D > 3), these equations are given by 

d 



dr 



r n — n + n 

dr 

d 



7T, 



a/3 



71 



ft i' 



2 V pvnP d a u p 



(4a) 
(4b) 



possible other second-order terms 24| can be neglected for the purpose of a linear stability 
analysis. In Eqs. @, the comoving derivative is denoted by u^d^ = d/dr. The relaxation 



times for the bulk viscous pressure and the shear stress tensor are denoted by rn and r n , 
respectively. The coefficients (, t] are the bulk and shear viscosities, respectively. We also 
introduced the symmetric rank-four projection operator 

1 1 

The shear stress tensor is traceless 7r^ = and orthogonal to the fluid velocity it /x 7T Ml/ = 0. 
The stability and causality of a relativistic dissipative fluid with bulk viscous pressure 

Thus, for the sake of simplicity, we shall first ignore 
the effects from bulk viscous pressure and discuss the properties of the fluid-dynamical 
equations of motion including only shear viscosity. The interplay between shear and bulk 
viscosity will be discussed afterwards. 



only have been investigated in Ref. 28 



A. Shear viscosity only 

For convenience, we introduce the following parameterization: 

rj = as , (6a) 

TV = _5_6 = ^, (6b) 
e+P T K ' 

where s and T are the entropy density and the temperature, respectively. From the second 
equation we obtain r n (e + P)/rj = b. The parametrization fl6]) is motivated by the leading- 
order results for the causal shear viscosity coefficient and the relaxation time obtained in 



Ref. (3l| where the relation tv — r]/P was found. For a massless ideal gas equation of state, 
e = (D — 1)P, this result is reproduced by choosing b = D. 

In this section, we discuss the stability of second-order relativistic fluid dynamics in the 

let us introduce a perturbation ~ e ^t-ikx arounc [ 
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local rest frame. Following Ref. 
the hydrostatic equilibrium state, 

e = e + 5e e wt ~ lkx , (7a) 
vr^ = Trft" + e iult - ikx , (7b) 
= u& + 5u» e iLut - ikx , (7c) 

where Eq = const., Hq U = 0, and Uq = (1, 0, 0, . . .), respectively. In the linear approximation, 
the velocity perturbation has no zeroth component, 

Su" = (0, Su 1 , 5u 2 , . . . , Su - 1 ) , (8) 
6 



because u^u^ = 1. Moreover, in the local rest frame, 5n = on account of the orthogonal- 
ity condition u^ir^ = 0. Since 7r M!y is traceless, Sti^ " 1 ^ ^ is not an independent variable. 
Taking all of this into account, the linearized fluid-dynamical equations can be written as 



AX = 



(9) 



where 



X = {5e,5u\5^\5u\5n 12 ,...,5u D -\5Ti l{D - l \ 
5vr 22 , 5tt 33 , . . . , 6^ D - 2 ^ D ~ 2 \ «5vr 23 , 5vr 24 , . . . , 



5n^ D - 1 \57r 3 \...J^ D -^ D ' 1 Y 



The matrix A is expressed as 



with 
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/TO \ 

B 

G C 

\ E J 



T 



G 

C 
E 



iu fx ^ 
—ike 2 fz —ik 

o r / y 



£ = diag(5 ,...,B )(D-2) 



x(D-2) , -B 




r 2 



o r 2 o 



(D-3)x3 
diag(J, ...,/) (D-3)x(D-3) , 

diag(/, . . . , /)i(o_ 2 )(D-3)x|(D-2)(D-3) > 



(10) 



(11a) 

(lib) 

(11c) 

(lid) 

(lie) 



where c s = y/dPjde is the velocity of sound. Here, we introduced the abbreviations 



h = iu(e + P) 
r x = — ikrj , 



f 1 = -,k(c + P), 



D — 1 
2 

»7 • 



L> - 1 



For nontrivial solutions of Eq. fl9]), the determinant of the matrix A should vanish. This 
leads to the following conditions for the dispersion relations u(k): 

f = , (12a) 



detB = (detB ) D ~ 2 
det T = det -%k c 2 s f 2 —%k 

{ o r / j 

Equation (I12a|) gives a purely imaginary frequency 







(12b) 
(12c) 



CO 



(13) 



which corresponds to a nonpropagating mode. The degeneracy of this mode is (D — 3)[1 + 
(D-2)/2\. 

Equation (112bj) leads to a complex frequency, 



2r^ 



corresponding to two propagating modes, if k is larger than the critical wavenumber 

k c = 



P Vb 



(14) 



(15) 



4 1] r n 2r n 

Following Ref. [36|, we shall call these modes shear modes. There are in total 2(D — 2) shear 
modes. 

Equation (|12c[) gives the same dispersion relation as Eq. (16) of Ref. [28|, after replacing 
2(D — 2)rj/(D — 1) with £ . Introducing the sound attenuation length in D space-time 
dimensions 

_ 2{D-2) V 2(D - 2) r n 
s ~ D-l e+P D-l b ' 
the analytic solution in the limit of small wavenumber k is 

i 



(16) 



UJ 



±kc s + i-^-k 2 



(17) 



while for large wavenumber we obtain 



i 



~1 -1 



1 + 



W = < 



±fc c s4 /l + At + 



7~ttC 2 2T-k 



T C 2 



r. 



(18) 
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k/T k/T 



FIG. 1: The real parts (left panel) and the imaginary parts (right panel) of the dispersion relations 
for the sound modes (full lines) and the nonpropagating mode (dashed line) obtained from Eq. 
(|12c[) . The parameters are a = , b = 6, c 2 = 3 for the 3+1-dimensional case, D = 4. 

This corresponds to another nonpropagating mode and two propagating modes which we 
call sound modes in accordance with Ref. 36|]. All imaginary parts are positive and therefore 



the nonpropagating, as well as the shear and sound modes are stable around the hydrostatic 



equilibrium state. This fact is already known from the study of 

In order to discuss the issue of causality, we follow Ref. 
velocity defined as 



iiscock and Lindblom 



28 



38] 



and study the group 



dReu 

B ' = "¥- (19) 

For the two nonpropagating modes, Re a; = 0. Consequently, in order to discuss causality, 
we have to consider the behavior of the imaginary part 



28 



. Let us digress for the moment 



and consider the diffusion equation with diffusion constant Dq. There is a nonpropagating 
mode with dispersion relation u = iD k 2 . Moreover, it is known that the diffusion equation 
is acausal. Therefore, we conjecture that a k 2 dependence of any nonpropagating mode 
can be considered a sign of acausality. In our case, the nonpropagating modes are either 
independent of k, or have a weak k dependence (cf. Fig. [I]). According to our conjecture, 
we conclude that the nonpropagating modes do not violate causality. 

The dispersion relations resulting from Eq. ( I12cj) are shown in Fig. [TJ and the correspond- 
ing group velocity resulting from Eq. ( fl~9j) in FigfSJ The group velocity has a maximum for 
a finite value of k/T and approaches its asymptotic value {k — > 00) from above. For small 
values of b, it may thus happen that the group velocity becomes superluminal. Nevertheless, 
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FIG. 2: The group velocity {22]) for a = l/(4vr), D 
(dashed line), as well as b = 1.5 (dotted line). 



4, (? s = ±, and b = 6 (full line), 6 = 2 



in Sec. IIII CI we shall show that only the asymptotic value determines whether the theory as 
a whole is causal or not. The asymptotic value of the group velocity is 

dReu 



lim 




(20) 



u q, sound a/ 
a ' fc-s-oo OK 

Consequently, for the asymptotic group velocity of sound waves to be less than the speed of 



light, r n and T s should satisfy the following, so-called asymptotic causality condition: 



< 1 



1 



V 



< 



D 



(21) 



r„ s b r n (e + P) ~ 2(D-2) K 

This is similar to the causality condition for the group velocity in the case of bulk viscosity, 
Eq. (21) of Ref. j^sj. For conformal fluids, where c 2 s = 1/(D — 1), the condition (121]) simplifies 
to r s < (D — 2)t„./ (D — 1) or, equivalently, b > 2. For example, for the values of t] and r n 



deduced from the AdS/CFT correspondence |36|, |40j |4lfl , 77 = s/(4vr), r ff = (2 - ln2)/(27rT), 
the condition (12 ip is always satisfied because b = 2(2 — In 2) ~ 2.614 > 2. 

The dispersion relations for the shear modes resulting from Eq. ( 112bl) change their be- 
havior from nonpropagating to propagating at the critical wavenumber ([To']) , as shown in 
Fig. [3] It should be noted that a similar behavior is observed in the case of bulk viscosity, 
cf. Fig. 1 in Ref. 28(. For wavenumbers larger than k c , the (modulus of the) group velocity 
of the propagating mode is 



v 



k/k c 



g, shear 



where 



g, shear 



V 



(22) 



(23) 
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k/T k/T 



FIG. 3: The real parts (left panel) and the imaginary parts (right panel) of the dispersion relations 
for the shear modes obtained from Eq. (|12bp . The parameters are a = , 6 = 6, c;? = ^ for the 
3+1-dimensional case, D = 4. 

V 

4^ 
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FIG. 4: The group velocity (ESJ) for D = 4 , 6 = 6, = ±, and a = l/(4vr) (full line), a = 1/4 
(dashed line), as well as a = 1 (dotted line). 

is the asymptotic value of v g in the large- wavenumber limit. If the asymptotic causality 
condition ([21]) is satisfied, v™ BheBX < ^/(D - 1)(1 - c 2 s )/2(D - 2). This is smaller than 1 for 
any value of c s and D > 3. However, near the critical wavenumber k c the group velocity 
diverges, as shown in Fig. HI From the definitions of k c , Eq. ( |T5|) . and the parameters a,b, 
Eqs. (EJ), we observe that k c /T = (2a\/b)~ 1 . The 1/a-scaling of k c /T for fixed 6 can be nicely 
observed in Fig. HI 

In Sec. IIII CI we shall show that the apparent violation of causality of the group velocity 
does not cause the theory as a whole to become acausal. The important issue is whether 




' k/T 

12 3 4 5 6 



11 



5 10 15 20 5 10 15 20 

k/T k/T 



FIG. 5: The real parts (left panel) and the imaginary parts (right panel) of the dispersion relations 
for the sound modes obtained from Eq. (|12cj) . The parameters are a = , 6=1, (? s = \ for the 
3+1-dimensional case, D = 4. 




k/T k/T 

FIG. 6: The real parts (left panel) and the imaginary parts (right panel) of the dispersion relations 
for the shear modes obtained from Eq. (|12bj) . The parameters are a = , 6 = 1, c 2 s = ^ for the 
3+1-dimensional case, D = 4. 

the asymptotic causality condition is fulfilled. If yes, the theory is causal. 

We remark that, in the local rest frame, the stability of the system of fluid- dynamical 
equations is not affected if we choose a parameter set which violates the asymptotic causality 
condition §ZT§, for instance a conformal fluid in D = 4 dimensions and 6=1. This is 
demonstrated for the sound modes in Fig. [51 and for the shear modes in Fig. El 
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B. Competition of bulk and shear 



The question we would like to answer in this section is whether the problem of the 
divergent group velocity can be removed by adding bulk viscosity to the discussion. For the 
sake of simplicity, we consider only the 2+1-dimensional case, i.e., D = 3. Similarly to Eqs. 
OH]), we introduce the parametrization 

c 



( = ais 



e + P 



(24) 



As before, the equations of motion fl4]) have to be linearized, yielding Eq. ([9]), where now 

X = (fe, 5u x , 8ir xx , 8u y , 5n x \ 5U) T , (25) 



and 



.4 



/ iu -ik{e + P) 

— ik c 2 iu (e + P) —ik 

— ik 7] iu r n + 1 

iu{e + P) 



\ -ikC, 







— ik 

ikt] iujr^ + 1 






—ik 




iu t u + 1 ) 



(26) 



Then, the dispersion relations are given by solving the following equations: 



k 2 r] + iu(l + iur n )(e + P) = 0, (27a) 



iu 



k 2 (l + iuT n )r] + (1 +iur^) [iuk 2 ( + (1 + iu r n )(e + P)(c 2 s k 2 



u 



. (27b) 



The dispersion relations resulting from sound and bulk viscous modes, Eq. (127b|) . are 

T 



u = < 



± [Aaaic 2 s (b + b x + bhc 2 ) - (a + at + abc 2 s + a^c 2 ,) 2 } 1 ^ j 

b 



(28) 



iT 



'1 

1 b h s 2(6 + h + bb lC 2 ) V«i&i ' ab 



for large k, and 



u 



i 



(29) 



±c 2 k, 
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for small k. 

Thus the asymptotic causality condition reads 

On the other hand, the equation for the shear modes, Eq. (127a[) . is the same as Eq. ( 11 2b>|) 



and hence the corresponding group velocity again shows a divergence. Thus, the inclusion 
of bulk viscosity does not solve the problem of the divergent group velocity. 

III. STABILITY IN LORENTZ-BOOSTED FRAME 



The discussion of causa 
boosted frame in Ref. 



28 



ity and stability in the case of nonzero bulk viscosity in a Lorentz- 
] has shown that causality and stability are intimately related. 
Relativistic dissipative fluid dynamics becomes unstable if the group velocity exceeds the 
speed of light. If this is still true in the case of nonzero shear viscosity, the divergence of the 
group velocity found in the rest frame may induce an instability in a moving frame. In order 
to investigate this question, we consider the stability of the hydrostatic state observed from 



28j. In this section, we restrict our investigations 



a Lorentz-boosted frame, following Ref. 
to the case D = 4. 

We consider a frame moving with a velocity V with respect to the hydrostatic state. 
Then, the total fluid velocity u' M is given by 

«"■=! lv Vlv * |< (3i) 

Vlvn 7vP|| + Q_ 



where 7y = — V 2 , Pn = nn T , and Q± = 1 — P\\, with n = V/\V\. We consider the two 

cases where the direction of the Lorentz boost is parallel and where it is perpendicular to 
the direction of propagation of the perturbation; the latter we take to be the x direction. 

A. Boost along the x direction 

The perturbation of the fluid velocity is given by 

u'» = u'^ + 5u'^ e luJt ~ ikx , (32) 
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where 



u'f = lv (1,V,0,0), 
Su'* 1 = (V^vSu x ,j v Su x ,6u y ,6u z 



(33a) 
(33b) 



where <5w M is the velocity perturbation in the local rest frame. The linearized fluid-dynamical 
equations are again given by Eq. (jUJ), with 



X = (8e, 5u x , 5tt xx , 5u\ 5n x \ 8u z , 5tt xz , 5tt v \ 5ix yz ) T 



(34) 



and 



.4 



/ Ti \ 

B x 

G x d 

\ E x ) 

The submatrices are given by 

( iuj{\ + V 2 c 2 s ) - ikV{\ + c 2 ) i[2uV - k(l + V 2 )} (e + P) i lv 2 V(uV - k) ^ 
iuV(l + c 2 ) - ik{V 2 + c 2 ) i [oj(1 + V 2 ) - 2kV] (e + P) i-fy 2 (coV-k) 

lir) lv (uV - k) rfF j 



(35) 



Tl = ly 



\ 



Bi = diag(B i,-B i) , B, 



01 



ijv(u- kV)(e + P) i{uV-k) 
irj"/y((juV — k) F 



Gx = ^0 
Ci = Ei — F 



-irj^viuV — k) 



Here we abbreviated 



Obviously, 



F = i-f V (u - kV)r n + 1 



detA = detTi x detBi x F 2 



From F 2 = 0, we only obtain two trivial propagating modes 
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+ kV . 



(36a) 
(36b) 

(36c) 
(36d) 

(36e) 
(37) 

(38) 



The group velocity is v g = V, which implies that these modes correspond to the nonpropa- 
gating modes in the LRF. 
From det£>x = 0, we obtain 



[iT + ab-f V (kV - u)]{W - to) + aj v (kV - tu) 2 T = 



(39) 



corresponding to the shear modes. There are in total four modes satisfying this relation. 
The solutions are given by 
1 



i T - 2a(l - b)W-iv ± J-T 2 + AiakTV-fy 1 + Aa 2 bk 2 -fy 2 



2a(b - V 2 h v 
On the other hand, the sound modes result from 

c 2 s (e + P) \ 1 - ijvr v (kV - u) 



(40) 



V 2 + (V - l) 2 V lv + 1 



2kVu 



(V-l)V lv -l + V 2 u 2 -c-\u-kV) 



kV 



<l V V(l - V) 



. 



(41) 



In Fig. [71 the dependence of the group velocity on the wavenumber is shown for various 
values of the boost velocity V. The left panel shows the behavior of one of the shear modes 
and the right panel one of the sound modes. The parameter set used here is a = b — 
6, c 2 s — g, which satisfies the asymptotic causality condition. We observe that the divergence 
of the group velocity of the shear mode in the rest frame is tempered by the Lorentz boost 
to result in a peak of finite height. However, the group velocity may still exceed the speed of 
light in a certain range of wavenumbers. As we increase the boost velocity, the peak height 
diminishes, until the group velocity remains below the speed of light for all wavenumbers. 
However, further increasing the boost velocity leads to an acausal group velocity in the 
sound mode. 

Although the group velocity of the shear or the sound mode may exceed the speed of 
light, as long as the asymptotic causality condition is fulfilled, the theory is still stable. This 
is demonstrated in the left panel of Fig. [HI where the imaginary parts of the modes are 
shown for the parameter set a = b = 6, c 2 — |. We observe that all imaginary parts are 
positive, indicating the stability of the theory. 

In contrast to the rest frame, where the theory is stable even for parameters which violate 
the asymptotic causality condition (|2T|) . this is no longer the case in a Lorentz-boosted 
frame. In the right panel of Fig. [HJ the imaginary parts of the modes are calculated with 
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FIG. 7: The group velocity calculated for one of the shear modes (left panel) and one of the sound 
modes (right panel). We set a = l/(4-7r),& = 6, (? s = 1/3. The solid line is for a boost velocity 
V = 0.05, the dashed line for V = 0.4 and the dotted line for V = 0.99, respectively. 
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FIG. 8: The imaginary parts of the dispersion relations for a boost in x direction with velocity 
V = 0.9. The left panel shows the results for the parameter set a = b = 6, c^ = |, which 
fulfills the asymptotic causality condition, while the right panel is for a = b = 1, c^ = g, which 
violates this condition. The dashed lines are for the shear modes, while the solid lines are for the 
sound modes. 



the parameter set a = j-, b = 1, c s = |. Now one observes the appearance of negative 
imaginary parts, indicating that the theory becomes unstable. 
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B. Boost along the y direction 



Now we consider a Lorentz boost along the y direction. The perturbation of the fluid 
velocity is given by 

(42) 



where 



n'" = u'f + Su'" e luJt ~ ikx 



u'f = 7^(1,0,^,0), 
8u ,fi = (V lv 5u\5u x , lv 5u y ,5u z ) . 



(43a) 
(43b) 



Similarly to the preceding discussion, the linearized fluid-dynamical equations take the form 
(jHJ), where the matrix A is 

/ T 2 H x H 2 \ 
H 3 B 2 H 4 H 5 
G 2 Hq C 2 
\ H 7 E 2 J 



A 



(44) 
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with 



To 



Hi 



Ho 



Bo 



#4 
G 2 
H 7 



( tU lv {l + C 2 g V 2 ) _ iklv ( £ + p) N 

—ike 2 iwyv(£ + P) —ik 

y -^ikr) F 1 

i 2iuV(e + P) 7 j> -ikV ^ 







T 



\ 

ituV 2 




V o 

/ iu^ v (l + V 2 )(e + P) -ik 



iuV 




-ikV-f V (e + P) o\ 
iuV^yfj 





V 



-ikjyT] 






Fi 





0/ 





\ 





iuj^y(e + P) —ik 
-ikr) Fi ) 



iuV 
l%k^ v r] 
zW7^ 



H 5 
He- 

C 2 : 



iuV 



T 



\iuV^\r] 



E 2 = F\ . 



(45a) 



(45b) 



(45c) 



(45d) 



(45e) 

(45f) 
(45g) 
(45h) 



Here we abbreviated 

F l = iuj^ivT^ + 1 . 

The condition detA = leads again to the following nine modes: three nonpropagating 
modes, four shear modes and two sound modes. 

The nonpropagating mode has almost the same form as that in the LRF, 



CO 



(46) 



The shear modes are given by the solution of the following equation 



k 2 r] + [V 2 7y?7w + (e + P)(i - jyT^u)] = , 



(47) 
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FIG. 9: The real and imaginary parts for the dispersion relations of the shear modes (dashed lines) 
and sound modes (solid lines), for a Lorentz boost in y direction. We use a = b = 6, c 2 = 
3, V = 0.9 in the 3+1-dimensional case. 



and the solutions are given by 

1 



2a(b - V 2 ) lv 



T ± V-T 2 + Aa 2 bk 2 - 4a 2 PV 2 



(48) 



We find that the critical wavenumber is now given by k c = T / (2a\Jb — V 2 ), below which the 
shear modes become nonpropagating modes. 

On the other hand, the sound modes and another nonpropagating mode result from 

3c 2 s (s + P)(-i + lv r^){k 2 + V^W) 
+ -y v u {Ak 2 r] + lv u [3i(e + P) + W 2 -f V r]u - 3{e + Phv^w] } = . (49) 

The real and imaginary parts of this dispersion relation are calculated with a parameter set 
satisfying the asymptotic causality condition. The results are shown in Fig. [9j One observes 
that the real parts are symmetric around u = 0. This symmetry is due to the fact that the 
direction of the Lorentz boost is orthogonal to the direction of the perturbation. The critical 
wave number k c where the shear mode changes from nonpropagating to propagating mode 
can be clearly seen. The imaginary parts are seen to be positive. We confirmed that the 
imaginary parts become negative if we use a parameter set which violates the asymptotic 
causality condition. 

C. Causality of wave propagation 

In the preceding discussion we have seen that the theory is stable if the asymptotic 
causality condition is fulfilled. The reverse is in general not true, as the discussion in the 
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local rest frame has shown, since a stable theory may also violate the asymptotic causality 
condition. However, the discussion in the Lorentz-boosted frame has revealed that the 
stability of a theory is contingent upon whether the asymptotic causality condition is fulfilled. 

In this section, we shall show that the causality of the theory as a whole is guaranteed if 
the asymptotic stability condition is fulfilled. The group velocity may become superluminal, 
or even diverge, as long as this apparent violation of causality is restricted to a finite range 
of momenta. The argument leading to this conclusion is analogous to that of Sommerfeld 



and Brillouin in classical electrodynamics 42|, |43[. For instance, in the case of anomalous 



dispersion the group velocity may become superluminal, but the causality of the theory as 
a whole is not affected. 

The change in a fluid-dynamical variable induced by a general perturbation is given by 

SX(x,t) =X) / duiXj^e**-*^* , (50) 
i 

where 8X(x, t) stands for Se, <5w M , and Stt^. The index j denotes the different modes, i.e., 
the shear modes, the sound modes etc. The function kj(u) is the inverted dispersion relation 

ujj(k) of the respective mode. The Fourier components are given by 

. , i r°° 

Y^SXj(u) = — / dt6X(0,t)e- ilJt . (51) 

We assume that the incident wave has a well-defined front that reaches x = not before 
t = 0. Thus 5X(0,t) = for t < 0. This condition on SX(0,t) ensures that ^2jSXj(u) is 
analytic in the lower half of the complex u plane 42j. On the other hand, in Sec. Ill Al we 
have found that the group velocity of the shear modes diverges for certain values of k. These 
divergences correspond to singularities in the complex uj plane. However, if the asymptotic 
causality condition is fulfilled, the imaginary part of the dispersion relation is always positive, 
i.e., the singularities only appear in the upper half of the complex u plane. In this case, the 
system is also stable. On the other hand, if the asymptotic causality condition is violated, 
the singularities may appear also in the lower half-plane, i.e., for negative imaginary part of 
the dispersion relation, and the system is unstable. 

We shall now demonstrate that the divergences in the group velocity do not violate causal- 
ity as long as the asymptotic causality condition is satisfied, i.e., as long as the asymptotic 
group velocity remains subluminal. To this end, we compute Eq. ( 150]) by contour integration 
in the complex u plane. To close the contour, we have to know the asymptotic behavior of 
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the dispersion relations. In our calculation, we found that the real part of the dispersion 
relation at large k is proportional to k [see Eq. (fT8]) ]. with a coefficient which is the large-/c 
limit of the group velocity, i.e., v^, 

lim Re Uj(k) = v^k . (52) 



k— >oo 



Then, in the large- A; limit, the exponential becomes 



exp[iut — ikj(ou)x] — > exp 



gj 



(53) 



In the case x > t, we have to close the integral contour in the lower half plane. If the 
asymptotic causality condition is fulfilled, there are no singularities in the lower half plane, 
and Eq. ( 1501) vanishes. On the other hand, the contour should be closed in the upper half 
plane if x < t. Then, because of the singularities, Eq. ( 150]) may have a nonzero value. 
However, as long as we choose a parameter set for which the asymptotic group velocity is 
smaller than the speed of light, i.e., for which the asymptotic causality condition is fulfilled, 
the signal propagation does not violate causality, since the locations x where the disturbance 
has travelled lie within the cone given by v^j which, in turn, lies within the lightcone, q.e.d. 

To conclude this section, we have shown that the asymptotic causality condition not only 
implies stability in a general (Lorentz-boosted) frame, but also causality of the theory as a 
whole. 

IV. CHARACTERISTIC VELOCITIES 

So far, we have analyzed the causality and stability of relativistic dissipative fluid dynam- 
ics with shear viscosity using a linear stability analysis. However, there is another possibility 
to analyze causality, namely by studying the characteristic velocities. For the sake of sim- 
plicity, we consider the 2+1-dimensional case with shear viscosity only. The fluid-dynamical 
equations can be written in the following form: 

{Ai b d t + A* ab d x + Al b d y ) Y b = B a , (54) 

where Y b T = (e, u x , u y , ir xx , ir xy ) and = (0, 0, 0, 7r xx , n xy ). The expressions for the 



components of A are given in the Appendix. Then, as discussed in Ref. 
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, the characteristic 
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velocities are defined as the roots of the following equations, 



det(^.A* - A x ) = , 
det(v4* - A v ) = . 



For the case of bulk viscosity, see Ref. 

For the sake of simplicity, we consider 
characteristic velocities are given by 



1, 0, 0) and ix 2 



71 



(55a) 
(55b) 

0. Then, the 





± 

± 



(56) 



Interestingly, the second velocity is identical to the asymptotic group velocity ( 1231) for the 
shear modes and the third velocity is the same as the asymptotic group velocity ([20]) for the 
sound modes (since D — 3). As a matter of fact, if the asymptotic causality condition ( F2~T|) 
is satisfied, the velocity (|56|) is smaller than the speed of light. 

In Fig. [TOl we show the b dependence of one of the five characteristic velocities. We set 
u M = (v^5/2, 1/2, 0), 7c xx = TT xy = 0, and c 2 s = 1/2. The velocity exhibits a divergence at 
small values of b, and thus exceeds the speed of light. This divergence occurs also for at 
least one other characteristic velocity. As far as we have checked numerically, in order to 
satisfy causality, one should use a value of b which is larger than about 2. This condition is 
consistent with the asymptotic causality condition (l2Tj) . 



V. CONCLUDING REMARKS 



In this work, we have discussed the stability and causality of relativistic dissipative fluid 
dynamics, based on a linear stability analysis around a hydrostatic state. Following the usual 
argument, we calculated the group velocity from the dispersion relation of the perturbation. 
We found that the group velocity diverges at a critical wavenumber k c . The appearance 
of the divergence is independent of the dimensionality of space-time and can be removed 
neither by tuning the parameters of the theory nor by adding bulk viscosity to the discussion. 

Nevertheless, in the rest frame of the background this acausal group velocity does not 
cause the fluid to become unstable. Moreover, investigating causality and stability in a 

23 



-J b 
5 



FIG. 10: One of the five characteristic velocities determined from the roots of Eqs. (155|) . The left 
panel is for v x and the right panel is for v y . We set = (y/b/2, 1/2, 0), ir xx = ir xy = 0, and 
c\ = 1/2. 

Lorentz-boosted frame, we found that the fluid-dynamical equations of motion are stable, if 
we choose parameters which satisfy a so-called asymptotic causality condition. They become 
unstable if this condition is violated. In this sense, the problems of acausality and instability 
are still correlated even in the case of shear viscosity, as was already found for the case of 
bulk viscosity 



We have then demonstrated that the causality of the theory as a whole is guaranteed 
if the asymptotic causality condition is fulfilled. Therefore, a superluminal group velocity 
in a finite range of momenta can cause the theory neither to become acausal nor unstable. 
Finally, we studied the characteristic velocities and found a violation of causality for small 
values of r n (e + P)/r), but not for values which satisfy the asymptotic causality condition. 

The asymptotic causality condition requires that the ratio T n /r s is sufficiently large, i.e., 
that the time scale r w over which the shear viscous pressure relaxes towards its NS value is 
not too small compared to the sound attenuation length T s ~ f]/{ £ + P) — v/{Ts)- This is 
an important finding for practitioners of fluid dynamics, who frequently consider tv and the 
shear viscosity-to-entropy density ratio 7]/s to be independent from each other. We have 
demonstrated that this is not the case if one wants the theory to remain causal. Therefore, 
second-order theories of relativistic dissipative fluid dynamics are not automatically causal 
by construction. Our findings also illuminate why NS theory violates causality from a 
different perspective, because there T n — > while 77 remains non-zero. 
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Appendix A: Matrix elements in Eq. (I54|) 

The fluid-dynamical equations can be expressed in the form ( 154)) . Let us parameterize 
the velocity of the fluid as = (cosh 9, sinh 9 cos 0, sinh# sin 0). The matrix elements of 

Kb are 

A\i — ( c s + l) s i nn $ cosn # cos , 

A x 12 = ^sech 3 6> | 2 sinh 2 9 [(2w + ir xx ) sin 2 + 3w cos 2 - ix xy sin cos 0] 

_ V 

+ 7/7sinh 4 #(cos(20) + 3) + w + tx xx } , 

A x 3 = sech 3 j sinh 2 6>cos0 [(w - n xx ) sin (ft + n xy cos0] + w; sinh 4 # sin cos + ir xy \ , 

A X A = tanh 9 cos , 

A x 5 = tanh 9 sin , 

A x 21 = (c 2 + l) sinh 2 9 cos 2 + c 2 , 

A 22 = 2u> sinh 9 cos , 

^24 — ^35 — 1 > 

^31 = ( c s + l) sin h 2 # sin cos , 

A x 2 = w sinh 9 sin , 

A33 = w sinh # cos , 

= sech 2 # I sinh 4 9 cos 2 [77 + T n n xx cos(20) - t^tx xx + tvtt 3 * sin(20)] 

+ sinh 2 9 [2(7/ - TV 71-^) cos 2 + 77 sin 2 0] + 77 j , 

A43 = -27V tanh 2 cos 2 [sinh 2 9 cos 0(tt C!/ cos - 7T m ' sin 0) + n xy ] , 

^44 = A x 5 = tv sinh 6 1 cos , 
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= tanh 2 cos r _ 2 sinh 2 gin 3 + ^ gin cog 2 + ^ CQS 

2(sinh 9 cos 2 + 1) 

+ sinh 4 9 sin 2 (20) (tt xs/ cos $ - 2ir xx sin 0) - 271-** sin - 271-^ cos } , 
= -sech 2 {2 sinh 4 9 cos 2 [77 - t^tt** cos(20) + 7-^ - r^sin^)] 
+ sinh 2 9 [( V + tvO cos(20) + 377 + T n Ti xx - T n 7T xy sin(20)] + 2^ } . 

The matrix elements of A l ah are given by 

4i = \ Ws + 1) cosh(2^) - c 2 + 1] , 



t 2 sinh 9 r 2 , 2 

A 9 = rr i sinh 6 l cos0 2wcos - 

(sinh 2 # cos 2 0+1) 1 
+ w sinh 4 9 cos 5 + (w + tt xx ) cos + ir xy sin 

A\ , = 2 sinh 6 1 w sin + 



7r xx sin 2 — 7T xy sin cos 



L 13 - * Olllll </ U/ Dill V, T . ,2/) T 

\ sinh t7 cos^ • 

/11 



, cos(20) 

14 ~~ v 2 



csch 9 + cos 2 1 
A t _ sin(20) 



csch 2 # + cos 2 ' 
A 2i = { c l + l) simi 61 cosn 9 cos 

^31 = ( C s + l) Sinn ^ COSn Sil1 

sech 3 # 



A\ 2 = — - — I 2 sinh 2 9 [(2w + ir xx ) sin 2 + 3w cos 2 - 7r X2/ sin cos < 

+ w sinh 4 [cos(20) + 3] + 2w + 2^ } , 
A 23 = sech 3 9 {sinh 2 9 cos [w sinh 2 9 sin + (w - ir xx ) sin + n xy cos 0] + tt xv } 



A 24 = tanh 6* cos , 

A 25 = tanh 9 sin , 

serin 3 r 

AL = ,\ sinh 2 fl \(w + ?,ti xx ) sin cos + 3vr^ sin 2 + 27^ cos 2 

(sinh 2 cos 2 0+l) 2 ^ L 

+ sinh 4 [3(w + tt**) sin cos 3 + (w + 50 sin 3 cos + 2^ sin 4 + tt X2/ cos 4 

+ — sinh 6 9 [10 sin(20) + sin(40)] \(w - ti xx ) cos(20) +w + tt xx - n xy sin(20)] 
16 

+ w sinh 8 sin cos 5 + 71^ 
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4s = 8 ( sinh 2g cos 2 0+1 ) { Sinh4 9 M W + C °< 2 ^ + (*" " "0 COS (^) + 21W 

- 9ir xx + 10vr xy sin(20) + ix xy sin(40)] + 4 sinh 2 [6w + 2vr xx cos(20) - 47T** 
+ 3vr xy sin(20)] - 4w sinh 6 cos 2 [cos(20) - 3] + 8w - 8ir xx } , 

tanh # sin (sinh 2 # sin 2 + l) 

34 sinh 2 6 1 cos 2 0+1 

, tanh 6> cos r . 2 . . /i 

2 smh cos z + 2 i> J 

A\ 2 = tanh 6 1 cos | sinh 2 6 1 {2 sin [(77 - r n 7i xx ) sin + t^tt^ cos 0] + 77 cos 2 0} + 77 - 2r 7r 7i 
A43 = - tanh 9 j sinh 2 6 1 cos 2 [(77 - 2r 7r 7r :i::!: ) sin + 2r 7r 7r x ' y cos 0] + 77 sin + 2r n 7i xy cos j 



^44 = ^55 = TV cosh 



A\ 2 = ( - 2 sinh2 9 l sin l~ 2 V + W x cos(20) + ?>t^ xx ] + 2r^7r^ cos 3 

4 sinh 6 cos 2 + 41 

+ sinh 4 9 sin 2 (20) [(77 - 2r 7T n xx ) sin + 2r 7r 7T xl/ cos 0] + 4(77 - t^ xx ) sin - 4r 7r 7r I ' s ' cos q 
A53 = tanh 9 j sinh 2 [77 cos 3 + Twr^ sin sin (20) - 2r 7r 7r :rs/ sin cos 2 0] 
+ (77 + T n 7r xx ) cos - T n ii xy sin j . 

The matrix elements of A v ab are 

A\ x = (c 2 + l) sinh 9 cosh 9 sin , 

^21 = ( c s + l) sm h 2 6* sin cos , 
sech 3 $ r 

^ 2 = ^ sinh 2 ^ [(7i; + 37r a;a; )sin0cos0 + 37r^sin 2 + 27r^cos 2 

(sinh 2 cos 2 + l) 2 I L 

+ sinh 4 [3(iu + ir xx ) sin cos 3 + (to + 57r xx ) sin 3 cos + 2^ sin 4 + w xy cos 4 
L sinh 6 0[1O sin(20) + sin(40)] [(to - tt xx ) cos(20) + w + n xx - n xy sin(20)] 
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+ w sinh 8 6 1 sin cos 5 + n xy j , 

= — = -{ sinh 4 0[4(w + 27r x:c )cos(20) + (7r :c:E -7i7)cos(40) + 2l7i; 

8 (sinh #cos 2 + 1) ^ 



A y 



97r xx + 10vr xy sin(20) + n xy sin(40)] + 4 sinh 2 9[6w + 2vr xa; cos(20) - 4vr a 
3tt^ sin(20)] - Aw sinh 6 cos 2 0[cos(20) - 3] + &w - 8n xx 
tanh 9 sin (sinh 2 9 sin 2 + l) 



sinh 2 9 cos 2 0+1 
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... tanh 9 cos r ,o nr , , s A 

AL = 5 <^ 2 - sinh 2 9 cos 20 - 3 \ 

15 2sinh 2 #cos 2 + 2 I J J 



A v 22 = w sinh 6* sin , 
A 23 = w sinh 9 cos , 

^25 = 1 ' 



A| x = (c 2 + 1) sinh 2 9 sin 2 + 

2 sinh [sinh 2 sin cos 0(7r 3:3: sin — ir xy cos 0) + ir xx cos + 7r xy sin 

^32 — 



(sinh 2 9 cos 2 0+1)' 



A33 = 2 sinh 6* ( w sin + 



sinh 9 cos 2 - 



sinh 2 6* sin 2 0+1 
j^y — 

sinh 2 9 cos 2 + 1 
sin(20) 



A v = - 

csch 2 6* + cos 2 ' 

A\ 2 = tanh 2 9 sin cos j sinh 2 9[2r] + cos(20) - r n ix xx + tvtt^ sin(20)] + 2rj - 2r n ir xx } 

A% = — j 2 sinh 4 cos 2 0[t? + r n 7T xx cos(20) - tvtt** + r^ xy sin(20)] 

+ sinh 2 0{r/[cos(20) + 3] + 2r^7r xy sin(20)} + 2r/} , 
^44 = ^55 = TV sinh 6* sin , 

A 52 = ^7— rr- D T——A sinh 2 0[(r 7r 7r^-r / )cos(40) + 9r / + 4r^'cos(20)-5r^ 

8 (sinh 9 cos z <p + 1) ^ 

- 87-^ sin cos 3 0] + 2 sinh 4 9 sin 2 (20) [77 + r^vr^ cos(20) - + sin(20)] 

+ 4[4?7 + t^tt** cos(20) - r^vr^ - r w vr XJ/ sin(20)] + 8r/csch 2 # | , 
A| 3 = tanh 2 sin [sinh 2 #sin(20) (71-^ sin - n xy cos 0) + ^ cos - 71-^ sin 0] , 

where we defined u; = £ + P. All other elements vanish. 
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